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Abstract. We investigate properties of Ehrhart polynomials for matroid poly- 
topes, independence matroid polytopes, and polymatroids. In the first half of 
the paper we prove that for fixed rank their Ehrhart polynomials are com- 
putable in polynomial time. The proof relies on the geometry of these poly- 
topes as well as a new refined analysis of the evaluation of Todd polynomials. 
In the second half we discuss two conjectures about the h* -vector and the co- 
efficients of Ehrhart polynomials of matroid polytopes; we provide theoretical 
and computational evidence for their validity. 



1. Introduction 

Recall that a matroid M is a finite collection of subsets T of [??] = {1, 2, . . . , n} 
called independent sets, such that the following properties are satisfied: (1) e JT, 
(2) if X e jc- and y C X then Y e (3) if U,V e T and \U\ = 1^^1 + 1 there 
exists X £ U\V such that VUx G In this paper we investigate convex polyhedra 
associated with matroids. 

One of the reasons matroids have become fundamental objects in pure and ap- 
plied combinatorics is their many equivalent axiomatizations. For instance, for a 
matroid M on n elements with independent sets J- the rank function is a func- 
tion ip: 2W -> Z where (p{A) := max{ |X| \ X C A, X e T}. Conversely a 
function tp: 2^"'^ Z is the rank function of a matroid on [n] if and only if the 
following are satisfied: (1) < ip{X) < |X|, (2) X C y =^ ip{X) < ip(Y), (3) 
ip{X UY) + ip{X n y) < ip{X) + ^(y). similarly, recall that a matroid M can 
be defined by its bases, which are the inclusion-maximal independent sets. The 
bases of a matroid M can be recovered by its rank function Lp. For the reader we 
recommend [31] or [40] for excellent introductions to the theory of matroids. 

Now we introduce the main object of this paper. Let B be the set of bases 
of a matroid M . li B = {ci, . . . , <Jr} G B, we define the incidence vector of B as 
bb '■= ^o^ii where is the standard elementary jth vector in M". The matroid 

polytope of M is defined as V{M) := conv{ \ B & B}, where conv(-) denotes the 
convex hull. This is different from the well-known independence matroid polytope, 
V'^{M) := conv{ ej \ I C B € B}, the convex hull of the incidence vectors of all the 
independent sets. We can see that P{M) C r^{M) and P(M) is a face of V^{M) 
lying in the hyperplane '^^^iXi = rank(M), where rank(Af) is the cardinality of 
any basis of M. 

Polymatroids are closely related to matroid polytopes and independence matroid 
polytopes. We first recall some basic definitions (see [40]). A function -0 : 2["1 — > R 
is submodular if '(l;{XnY)+ip{XUY) < ipiX) +i>(Y) foT all X,Y C [n]. A function 
is non- decreasing if tpi^) ^ V'(^) for all X C y C [n]. We say ip 
is a polymatroid rank function if it is submodular, non-decreasing, and V'(0) = 0- 
For example, the rank function of a matroid is a polymatroid rank function. The 
polymatroid determined by a polymatroid rank function ip is the convex polyhedron 
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(see Theorem 18.2.2 in [40]) in M" given by 

Independence matroid polytopes are a special class of polymatroids. Indeed, if 
is a rank function on some matroid M, then V^{M) = V{(p) [17]. Moreover, 
the matroid polytope V{M) is the face of Vlip) lying in the hyperplane X]"=i — 
ip{\n\). Matroid polytopes and polymatroids appear in combinatorial optimization 
[34], algebraic combinatorics [19], and algebraic geometry [21]. The main theme of 
this paper is the study of the volumes and Ehrhart functions of matroid polytopes, 
independence matroid polytopes, and polymatroids (from now on we often refer to 
all three families as matroid polytopes). 

To state our main results recall that given an integer fc > and a polytope 
-P C R" we define kV := {ka \ a. & V} (md the function i{V, k) := #(fc7' n Z"), 
where we define i{V, 0) 1. It is well known that for integral polytopes, as in the 
case of matroid polytopes, i{'P,k) is a polynomial, called the Ehrhart polynomial 
of v. Moreover the leading coefficient of the Ehrhart polynomial is the normalized 
volume of V, where a unit is the volume of the fundamental domain of the affine 
lattice spanned by V [36]. Our first theorem states: 

Theorem 1. Let r be a fixed integer. Then there exist algorithms whose input data 
consists of a number n and an evaluation oracle for 

(a) a rank function ip of a matroid M on n elements satisfying ip{A) < r for all A, 
or 

(b) an integral polymatroid rank function ip satisfying ip^A) < r for all A, 

which compute in time polynomial in n the Ehrhart polynomial (in particular, the 
volume) of the matroid polytope V{M), the independence matroid polytope V^{M), 
and the polymatroid 'P{ip), respectively. 

The computation of volumes is one of the most fundamental geometric opera- 
tions and it has been investigated by several authors from the algorithmic point of 
view. While there are a few cases for which the volume can be computed efficiently 
(e.g., for convex polytopes in fixed dimension), it has been proved that computing 
the volume of polytopes of varying dimension is #P-hard [8, 16, 25, 29]. Moreover 
it was proved that even approximating the volume is hard [18] . Clearly, computing 
Ehrhart polynomials is a harder problem still. To our knowledge two previously 
known families of varying-dimension polytopes for which there is efficient computa- 
tion of the volume are simplices or simple polytopes for which the number of vertices 
is polynomially bounded (this follows from Lawrence's volume formula [29]). Al- 
ready for simplices it is at least NP-hard to compute the whole list of coefficients of 
the Ehrhart polynomial, while recently [2] presented a polynomial time algorithm 
to compute any fixed number of the highest coefficients of the Ehrhart polynomial 
of a simplex of varying dimension. Theorem 1 provides another interesting family 
of varying dimension with volume and Ehrhart polynomial that can be computed 
efficiently. The proof of Theorem 1, presented in Section 2, relies on the geometry 
of tangent cones at vertices of our polytopes as well as a new refined analysis of 
the evaluation of Todd polynomials in the context of the computational theory of 
rational generating functions developed by [1, 2, 3, 4, 10, 13, 39, 41]. 

In the second part of the paper, developed in Section 3, we investigate algebraic 
properties of the Ehrhart functions of matroid polytopes: The Ehrhart series of 
a polytope V is the infinite series ^'^^Qi{'P,k)t''. We recall the following classic 
result about Ehrhart series (see e.g., [23, 36]). Let V C R" be an integral convex 
polytope of dimension d. Then it is known that its Ehrhart series is a rational 
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function of the form 



oo 



hi + hit + ■ ■ ■ + h*^_^t''-^ + h*^t' 



(1) 



The numerator is often called the h* -polynomial of V, and we define the coefficients 
of the polynomial in the numerator of Lemma 1, hQ + h^t + ■ ■ ■ + h'^^it'^^^ + h'^t'^, 
as the h*-vector of V, which we write as h*{V) := {hi, hi, /ij-n ^d)- 

A vector (cq, . . . , c^) is unimodal if there exists an index p, < p < d, such that 
Ci-i < Ci for i < p and Cj > Cj+i for j > p. Due to its algebraic implications, several 
authors have studied the unimodality of /i*-vectors (see [23] and [36] and references 
therein). It is well-known that if the Ehrhart ring of an integral polytope V, A{'P), 
is Gorenstein, then h*{'P) is unimodal, and symmetric [23, 36]. Nevertheless, the 
vector h*{V) can be unimodal even when the Ehrhart ring A{V) is not Gorenstein. 
For matroid polytopcs, their Ehrhart ring is indeed often not Gorenstein. For 
instance, De Negri and Hibi [15] prove explicitly when the Ehrhart ring of a uniform 
matroid polytope is Gorenstein or not. Two fascinating facts, uncovered through 
experimentation, arc that all /i*-vectors seen thus far are unimodal, even for the 
cases when their Ehrhart rings are not Gorenstein. In addition, when we computed 
the explicit Ehrhart polynomials of matroid polytopes we observe their coefficients 
are always positive. We conjecture: 

Conjecture 2. Let 'P{M) he the matroid polytope of a matroid M . 

(A) The h* -vector of'P{M) is unimodal. 

(B) The coefficients of the Ehrhart polynomial ofV{M) are positive. 

We have proved both parts of this conjecture in many instances. For example, 
using computers, we were able to verify Conjecture 2 for all uniform matroids up 
to 75 elements as well as for a wide variety of non-uniform matroids which are 
collected at [12]. We include here this information just for the 28 famous matroids 
presented in [31]. Results in [24], with some additional careful calculations, imply 
that Conjecture 2 is true for all rank 2 uniform matroids. Regarding part (A) of 
the conjecture we were also able to prove partial unimodality for uniform matroids 
of rank 3, meaning that the vector is non-decreasing up to a non-negative index, 
for large enough n. Concretely we obtain: 

Theorem 3. 

(1) Conjecture 2 is true for all uniform matroids up to 75 elements and all 
uniform matroids of rank 2. It is also true for all matroids listed in [12]. 

(2) Let V(U^'") be the matroid polytope of a uniform matroid of rank 3 on n 
elements, and let I be a non-negative integer. Then there exists n{I) G N 
such that for all n > n(I) the h* -vector of V^U^'"), (/ig, . . . ,h^), is non- 
decreasing from index to I . That is, h^ < hi < ■ ■ ■ < hj. 



Generating functions are crucial to proving our main results. Let V C M" be a 
rational polyhedron. The multivariate generating function of V is defined as the 
formal Laurent series in Z[[zi, . . . , 2„, z^^ , . . . , z^^]] 



where we use the multi-exponent notation z" = n"=i ^t'- ^ bounded, g-p 
is a Laurent polynomial, which we consider as a rational function g-p. If V is 
not bounded but is pointed (i.e., V docs not contain a straight line), there is a 
non-empty open subset [/ C C" such that the series converges absolutely and 
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uniformly on every compact subset of f7 to a rational function g-p (see [3] and 
references therein). If V contains a straight line, we set g-p = 0. The rational 
fimction gp E Q(^i, ■ • ■ , defined in this way is called the multivariate rational 
generating function of V D Z". Barvinok [1] proved that in polynomial time, when 
the dimension of a polytope is fixed, gp can be represented as a short sum of 
rational functions 

iei -1^5=1^ ^ 

where Si £ { — 1, 1}. Om' first contribution is to show that in the case of matroid 
polytopcs of fixed rank, this still holds even when their dimension grows. 

2.1. On the Tangent Cones of Matroid Polytopes. Our goal is to compute 
the multivariate generating function of matroid polytopes and independence ma- 
troid polytopes with fixed rank (later, in Section 2.2, we will deal with the case of 
polymatroids), and to do this we will use a crucial property of adjacent vertices. 
Let V be a vertex of V. Define the tangent cone or supporting cone of v to be 

Cp{v) := { V + w I V + £w G V for some e > } . 

To illustrate our techniques we will use a running example throughout this sec- 
tion. 

Example 4 (Matroid on K4). Let K4 be the complete graph on 4 vertices. Label 
the (2) = 6 edges with {1, . . . , 6}. Every graph induces a matroid on its edges where 
the bases are all spanning trees (spanning forests for disconnected graphs) [40]. Let 
M{K4) be the matroid on the elements {1, . . . , 6} with bases as all spanning trees 
of K4. The rank of M{K4) is the size of any spanning tree of K4, thus the rank 
oi M{K4) is 3. The 16 bases of M{K 4) are: {3,5,6}, {3,4,6}, {3,4,5}, {2,5,6}, 
{2,4,6}, {2,4,5}, {2,3,5}, {2,3,4}, {1,5,6}, {1,4,6}, {1,4,5}, {1,3,6}, {1,3,4}, 
{1,2,6}, {1,2,5}, {1,2,3}. 

Lemma 5 (See Theorem 4.1 in [21], Theorem 5.1 and Corollary 5.5 in [38]). Let 
M be a matroid. 

A) Two vertices es^ and eg^ are adjacent in ViM) if and only if Gbi ^ ^3-2 ~ 
ei — ej for some i,j. 

B) // two vertices ej^ and e/^ are adjacent in {M) then e/^ — e/^ G { e^ — 
Bj, Bi, — ej } for some i,j. Moreover if v is a vertex of V^{M) then all 
adjacent vertices of w can he computed in polynomial time in n, even if the 
matroid M is only presented by an evaluation oracle of its rank function ip. 

Let M be a matroid on n elements with fixed rank r. Then the number of 
vertices of V{M) is polynomial in n. We can see this since the number of vertices 
is equal to the number of bases of M, and the number of bases is bounded by ("), 
a polynomial of degree r in n. Clearly the number of vertices of V^{M) is also 
polynomial in n. It is also clear that, when the rank r is fixed, all vertices of either 
polytope can be enumerated in polynomial time in ??, even when the matroid is 
only presented by an evaluation oracle for its rank function ip. 

Throughout this section wc shall discuss polyhedral cones C with extremal rays 
{ri, . . . , r;} such that 

rfe G Ra-= { e; - e^, e^, -e^ \ i E[n], j E A} for k^l,...,l 

for some A C [n]. We will refer to Ra as the elementary set of A. Note that by 
Lemma 5 the rays of a tangent cone at a vertex e^ (corresponding to a set A C_ [n]) 
of a matroid polytope or an independence matroid polytope form an elementary set 
of A. Due to convexity and the assumption that are extremal, for each i G [n] 
and j E A at most two of the three vectors e^ — e^, e^, — e^ are extremal rays rfc of 



EHRHART POLYNOMIALS OF MATROID POLYTOPES AND POLYMATROIDS 



5 



C. This implies by construction, that considering all pairs — e-,- and or — e^, 
the number of generators of C is bounded by 

n\A\+n+\A\. (2) 

Recall a cone is simple if it is generated by linearly independent vectors and it 
is unimodular if its fundamental parallelepiped contains only from Z" [22]. A 
triangulation of C is unimodular if it is a polyhedral subdivision such that each 
sub-cone is unimodular. 




Figure 1. {2,3,5}, {2,3,4}, {1,3,6}, {1,3,4}, {1,2,6} and 
{1, 2, 5} are spanning trees of K4 that differ from {1, 2, 3} by adding 
one edge and removing one edge. 

Example 6 (Matroid on K4). The vertices e|2,3,5}, e|2,3,4}, e|i 3g}, e|i 34}, 
6{i,2,6} cind ^{1,2,5} are all adjacent to the vertex e^i 2,3}j see Figure 1. More- 
over, the tangent cone C-p(^M(K4)){^{i,2,3}) is generated by the differences of these 
vertices with e|x,2,3}- 

Cp(A/(K4))(e{l,2,3}) = ^{1,2,3} + cone {e{2,3,5} - e{i,2,3}:e{2,3,4} - e{i,2,3}: 

^{1,3,4} ~ 2,3}; ^{1^2,6} " 2,3} ; ^{1,2,5} " ^{1,2,3}} 

Lemma 7. Let C C R" be a cone generated by p extremal rays {ri, . . . , Vp} C i?^ 
where Ra is an elementary set of some A C [n] . Every triangulation of C is 
unimodular. 

Proof. Without loss of generality, we can assume {ri, . . . , r;} are generators of the 
form e^ — e-,- and {vi^i, . . . , r^} are generators of the form Bi or — for the cone C. 

It is easy to see that the matrix Tc := [ri, . . . , r;] is totally unimodular. Let Gc 
be a directed graph with vertex set [n] and an edge from vertex i to j if r^, = — ej 
is an extremal ray of C. We can sec that Gc is a subgraph of the complete directed 
graph Kn with two arcs between each pair of vertices; one for each direction. Since 
Tc := [ri, . . . , r;] is the incidence matrix of the graph Gc, it is totally unimodular 
[33, Ch. 19, Ex. 2], i.e., every subdeterminant is 0, 1 or -1 [33, Ch. 19, Thm. 9]. 
Therefore Tc [ri, . . . , r/, rj+i, . . . , r^] is totally unimodular since augmenting Tc 
by a vector e^ or — e^ preserves this subdeterminant property: for any submatrix 
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containing part of a vector e,; or —Bj perform the cofactor expansion down the 
vector ei or — when calculating the determinant. 

Since Tq is totally unimodular, each basis of Tq generates the entire integer 
lattice Z" nlin(C) and hence every simplicial cone of a triangulation has normalized 
volume 1. □ 

Lemma 8. Let C C R" be a cone generated by I extremal rays {ri, . . . , r;} C Rj^ 
where Ra is an elementary set of some A C [n], where dim(C) < n. The extremal 
rays {ri, . . . , r/} can be augmented by a vector r such that dim(cone{ri, . . . , r;, r}) = 
dim(C) + 1,, the vectors ri, . . . , r/, r are all extremal, and r e Ra- 

Proof. It follows from convexity that at most two of e; — ej, or — e-,- are extremal 
generators of C for z e [n] and j E A. There are at least n possible extremal 
ray generators, considering two of — e^, or — for each i S [n] and j G A. 
Moreover, all these pairs span M". Thus by the basis augmentation theorem of linear 
algebra, there exists a vector r such that dim(cone{ri, . . . , r;, r}) = dim(C) + 1 and 
ri , . . . , r; , r are all extremal. □ 

Lemma 9. Let r be a fixed integer, n be an integer, A C [n] with |A| < r and let 
C C M" be a cone generated by I extremal rays {ri, . . . , r;} C Ra where Ra is an 
elementary set of A. Then any triangulation o/conv({0, ri, . . . ,rj}) has at most a 
polynomial in n number of top- dimensional simplices. 

Proof. Assume dim(C) = n. Later, we will show how to remove this restriction. 
We can see that conv{0, ri, . . . , r/} C [—1, 1]"^ x A[„]\^ where 

[-1,1]-^ := {xeM-^ I |xj| < 1 je^jCM^ 

A[„]\^ := conv({ \ i e [n] \A}U {0}) C R["1\^. 

The volume of a d-simplcx convjvo, . . . , v^} is [22] 

(3) 

Thus the {n — |yl|)-volume of A[„]\^ is and the |A|-volume of [—1, 1]^ is 

2l^l. Therefore 

vol ([-1, 1]-^ X A[„]\^) = 2l^l (^^3^ = ^2^Mn - 1) • • • (n - |A| + 1). 

It is also a fact that any integral n-simplcx has n- volume bounded below by , using 
the simplex volume equation (3). Therefore any triangulation of conv{0, ri, . . . , r;} 
has at most 

2'^ln(n - 1) • • • (n - 1^1 + 1) < 2'^n(n - 1) • • • (n - r + 1) 

full-dimensional simplices, a polynomial function in n of degree r. 

Let dc := n — dim(C). If dim(C) < n, then by Lemma 8, {ri,...,ri} can 
be augmented with vectors {ri, . . . ,fdc} where S Ra for A above, such that 
dim(cone{ri, . . . ,ri,fi, . . . ,r(ij,}) = n and {ri, . . . , r;, ri, . . . , rd^, } are extremal. 
Moreover, 

dim(conv{0, ri, . . . ,r;}) < dim(conv{0, ri, . . . ,r;,fi}) < • • ■ 
< dim(conv{0, ri, . . . ,r;,ri, . . . ,rdc-i}) < dim(conv{0, ri, . . . ,r;,ri, . . . ,fdc}), 

that is, rfe ^ aff{0, ri, . . . ,r;,ri, . . . ,ffc-i} for 1 < fc < dc. 

Since Ffc ^ aff{0, ri, . . . , r;, ri, . . . , Vk-i}, any fuU-dimcnsional simplex in a trian- 
gulation of conv{0, ri, r;, , r^} must contain r^, see Figure 2. If not, then 
there exists a top-dimensional simplex using the points {0, ri, . . . , r;, ri, . . . , r^-i}. 
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Figure 2. If is not contained in the afRne span of 
{0, ri, Tp, , ffc-i} then every fuU-dimensional simplex 

must contain r^. 

but we know all these points lie in a subspace of one less dimension, a contra- 
diction. Therefore, a bound on the number of simplices in a triangulation of 
conv{0, ri, r;, , Vk} is a bound on that of conv{0, ri, . . . , r;, ri, . . . , Vk-i}- 

Thus, if dim(C) < n we can augment C by vectors ri, . . . ,fdc so that the cone 
C := conejri, . . . , rp, ri, . . . , r^c} is of dimension n and r;, e Ra for A above. 
We proved any triangulation of conv{0, ri, . . . , rp, ri, . . . , rrf^} has at most polyno- 
mially many full-dimensional n-simplices, which implies that any triangulation of 
convjO, ri, . . . ,rp} has at most polynomially many top-dimensional simplices due 
to the construction of the generators □ 

We have shown that for a cone C generated by an elementary set of extremal rays 
{ri, . . . , r;} C Ra for some AC [n], any triangulation of conv{0, ri, . . . , rj has at 
most polynomially many simplices. What we need next is an efficient method to 
compute some triangulation of conv{0, ri, . . . , r;}. We will show that the placing 
triangulation is a suitable candidate. 

Let V C M" be a polytope of dimension n and A be a facet of V and v e R" . 
There exists a unique hyperplane H containing A and P is contained in one of the 
closed sides of H, call it H~^. If v is contained in the interior of H~ , the other 
closed halfspace defined by H, then A is visible from v (see chapter 14.2 in [22]). 
The well-known placing triangulation is given by an algorithm where a point is 
added to an intermediate triangulation by determining which facets are visible to 
the new point [14, 22]. We recall now how to determine if a facet is visible to a 
vertex in polynomial time. 




Figure 3. Figure (a) shows A visible to v. Figure (b) shows A 
not visible to v. 



8 



JESUS A. DE LOERA, DAVID C. HAWS, AND MATTHIAS KOPPE 



Lemma 10. Let V C M" he a polytope given by t vertices {v^, . . . , v*} C R" and 
A C V be a facet of V given by q vertices {v^, . . . jV"?} C {v^, . . . ,v*}. If v ^ M" 
where w ^ V then deciding if A is visible to v can be done in polynomial time in 
the input {v^, . . . , v''}, {v^, . . . , v*} and v. 

Proof. Let z i Yl\=i t\\&t z G rclint(A). We consider the linear program: 







1 









.=1 ^=1 

< A < 1, Av+ (1 - A)z = x|. 

If (4) has a solution then there exists a point 5i E V between the facet A and v, 
hence A is not visible from v. If (4) does not have a solution, then there are no 
points of V between v and A, hence A is visible from v (see Lemma 4.2.1 in [14]). 
It is well known that a strict inequality, such as the one in (4), can be handled by 
an equivalent linear program which has only one additional variable. Determining 
if (4) has a solution can be done in polynomial time in the input [33]. □ 

Algorithm 11 (The Placing Triangulation [14, 22]). 

Input: A set of ordered points {vi, . . . , vt} G M". 
Output: A triangulation T of {vi, . . . , Vf } 



1 

2 
3 

4: 
5 
6 
7: 
8 
9 
10 
11 

12 
13 
14 
15 



^:={{vi}}. 

for each G {v2, . . . , vt} do 
Let B eT. 
Pi := {vi, . . . , Vi_i} 
if Vi ^ aff(Pi) then 
T' :=0 

for each G T do 
r := T' U U {vj}. 

else 

for each B E T and each (|-B| — l)-subset C oi B do 

Create and solve the linear program (4) with {Pi,C,Vi) to 
decide visibility of C to Vi . 
if C is visible to then 
r := T' U {C U {v,}} 
T := T' 
return T 



Indeed, Algorithm 11 returns a triangulation [22]. Wc will show that for certain 
input, it runs in polynomial time. Wc remark that there are exponentially, in n, 
many lower dimensional simplices in any given triangulation. But. it is important 
to note that only the highest dimensional simplices are listed in an intermediate 
triangulation (and thus the final triangulation) in the placing triangulation algo- 
rithm. 

Theorem 12. Let r be a fixed integer, n be an integer, A C [n] with \A\ < r, 
and let C C R" be a cone generated by extremal rays {ri, . . . ,r;} C i?^. Then the 
placing triangulation (Algorithm 11) with input {0, ri, . . . , r;} runs in polynomial 
time. 



Proof. By Equation (2) there is only a polynomial, in n, number of extremal rays 
{ri, . . . ,r;}. Thus, the for statement on line 2 repeats a polynomial number of 
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times. Step 5 can be done in polynomial time by solving the linear equation 
[ri, . . . ,ri_i]x = Vj. 

The for statement on line 7 repeats for every simplex D in the triangulation T, 
and the number of simplices in T is bounded by the number of simplices in the 
final triangulation. By Lemma 9 any triangulation of extremal cone generators in 
Ra with the origin will use at most polynomially many top-dimensional simplices. 
Hence the number of top-dimensional simplices of any partial triangulation T will 
be polynomially bounded since it is a subset of the final triangulation. 

The for statement on line 10 repeats for every simplex B and every {\B\ — 1)- 
simplex of B. As before, the number of simplices B is polynomially bounded, and 
there are at most n {\B\ — l)-simplices of B. Thus the for statement will repeat a 
polynomial number of times. 

Finally, by Lemma 10, determining if C is visible to can be done in polynomial 
time. Therefore Algorithm 11 runs in a polynomial time. □ 

Corollary 13. Let r be a fixed integer, n be an integer, A C [n] with \ A\ < r, and 
let C C R" be a cone generated by extremal rays {ri, . . . , r/} C R^. A triangulation 
of C can be computed in polynomial time in the input of the extremal ray generators 
{ri, . . . ,r;}. 

Proof. Let Vc conv{0, ri, . . . , r^}. We give an algorithm which produces a tri- 
angulation of Pc := conv{0, ri, . . . , rj} such that each full-dimensional simplex has 
as a vertex. Such a triangulation would extend to a triangulation of the cone C. 
This can be accomplished by applying two placing triangulations: one to triangu- 
late the boundary of Pc not incident to 0, and another to attach the triangulated 
boundary faces to 0. The algorithm goes as follows: 

1) Triangulate Pc using the placing triangulation algorithm. Call it 
T'. 

2) Triangulate Pc using the boundary faces of T' which do not con- 
tain V. 

Algorithm 14 (Triangulation joining to boundary faces). 

Input: A triangulation T' of Pc, given by its vertices. 
Output: A triangulation T of Pc such that every highest di- 
mension simplex of T is incident to 0. 
1: r := 

2: for each C where C is a — l)-simplex oi B E T' do 
3: if C is not a - l)-simplex oi A e T' where A ^ B 
then 

4: r :=ru{cu{o}} 

5: return T 




Figure 4. A triangulation T' of Pc can be used to extend to a 
triangulation T such that is incident to every highest dimensional 
simplex 
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By Theorem 12, triangulating Vc using Algorithm 11 can be done in polynomial 
time. Algorithm 14 indeed produces a triangulation of Vc- It covers Vc since every 
extremal ray generator of C is on some (dim(C) — l)-simplex. Moreover, T by 
construction has the property that the intersection of any two simplices of T is a 
simplex. Step 3 checks if C is on the boundary, since if C is on the boundary it 
will not be on the intersection of two higher-dimensional simplices. 

Step 2 repeats a polynomial number of times since any triangulation of Vc has 
at most a polynomial number of simplices, and each simplex B has at most n 
{\B\ — l)-simpliccs. Step 3 can be computed in polynomial time since again there 
are only polynomially many simplices B in the triangulation T' and at most n 
{\B\ — l)-simplices to check if they are equal to C. Hence, Algorithm 14 runs in 
polynomial time. □ 

Example 15. The tangent cone at the vertex := 2,3} on the polytope 
V{M{K4)) can be triangulated as: 

I {^{2,3,5} - es, e{2,3,4} - es, e{i_3^6} - ^s, e{i_3^4} - e^, e{i 2,6} - ^b] , 

{^{2,3,5} - ^B, e{i,3,6} - e_B, e{i 3 4} - es, e{i 2,6} - ^b, ^{1,2,5} - ^b} , 

{e{2,3,5} - es, e{2,3,4} - e^, e{i_3^4} - e^, e{i_2,6} - ^b, ^{1,2,5} - ^b] |- 

2.2. Polymatroids. We will show that certain lemmas from Subsection 2.1 also 
hold for certain polymatroids. Recall that the rank of the matroid M is the size of 
any bases of M which equals V'(W)- Our lemmas from Subsection 2.1 rely on the 
fact that M has fixed rank, that is, for some r S Z, r > 0, ^p{A) < r for all A CI [n]. 
We will show that a similar condition on a polymatroid rank function is sufficient 
for the lemmas of Subsection 2.1 to hold. 

Lemma 16. Let 2["1 — > N be an integral polymatroid rank function where 
'ip{A) < r for all A <Z [n], where r is a fixed integer. Then the number of vertices 
ofV{ip) is bounded by a polynomial in n of degree r. 

Proof. It is known that if ip is integral then all vertices of V{ip) are integral [40]. 
The number of vertices of V{ip) can be bounded by the number of non- negative 
integral solutions to xi + • • ■ + a;„ < r, which has ("^'^) solutions, a polynomial in 
n of degree r [37]. □ 

Lemma 17. Let ip: 2["1 — > N be an integral polymatroid rank function. If -v is 
a vertex of V{ip) then all adjacent vertices of v can be enumerated in polynomial 
time. Moreover if ijj(A) < r for all A C [n], where r is a fixed integer, then the 
vertices ofV{ip) can be enumerated in polynomial time. 

Proof. If V is a vertex of V{ip) then generating and listing all adjacent vertices to v 
can be done in polynomial time by Corollary 5.5 in [38]. If V'(^) ^ foi' a-U AC [n], 
where r is a fixed integer, then, by Lemma 16, there is a polynomial number of 
vertices for V{ip). We know that G M" is a vertex of any polymatroid. Therefore, 
beginning with 0, we can perform a breadth-first search, which is output-sensitive 
polynomial time, on the graph of 'P('0), enumerating all vertices of V{ip). □ 

What remains to be shown is that these polymatroids have cones like the ones in 
Subsection 2.1. We first recall some needed definitions from [38]. Let v, w £ R" and 
define A(v, w) := { ^ e [n] \ v, ^ w, } and cl(v) := { 5 | 5 C [n], Y.^es = ^i^) }. 
Let F = {/i, . . . , /|F|} be an ordered subset of [n] and Fi := {/i, . . . , fi}. li ip is a. 
polymatroid rank function then we construct v g R" where Vi = ip{Fi) — V'(-Fi-i) 
where vj = when j ^ F and one says F generates v. A classical result of Edmonds 
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[17] says that the set of vectors generated by all ordered subsets of [n] is exactly 
the set of vertices of V{ip). Now we can restate an important lemma. 

Lemma 18 (See Theorem 4.1 and Section 2 in [38].). Let tp be a polymatroid rank 
function. If v and w are vertices of the polymatroid V (tp) then either 

(i) |A(v, w)| = 1 or 

(ii) cl(v) = cl(w) and A(v,w) = {c,d} for some c,dG [n] where there exists 
some ordered set F ~ {/i, . . . , /|f|} which generates v with fk+i = d and 
fk ~ c for some integer k, 1 < k < \F\ — 1; moreover the ordered set 
F {/i, • ■ ■,fk-i,fk+i,fk,fk+2, ■ ■ • ,./|F|} generates w. 

Lemma 19. Let tJj be an integral polymatroid rank function andC the tangent cone 
of a vertex v of the polymatroid V^^), translated to the origin. Then C is generated 
by extremal ray generators {ri, . . . ,r/} C i?supp(v); where i?supp(v) an elementary 
set of supp(v). 

Proof. Let ip: 2l"l — > Z be a integral polymatroid rank function. Let v and w be 
adjacent vertices of the polymatroid 'P{ili). Using Lemma 18, if |A(v,w)| = 1 then 
w — V = Hbi where h is some integer and e^ is the standard ith elementary vector 
for some z S [n]. If /i < then certainly i £ supp(v), else i £ [n\. Thus w — v, a 
generator of C, is parallel to a vector in i?supp(v)- 

Let V and w be adjacent and satisfy (ii) in Lemma 18, where A(v,w) = {c, d}. 
Hence there exists an = {/i, . . . , f\F\} which generates v with fk+i ~ d and 
fk = c for some integer k, 1 < k < — 1; moreover the ordered set F := 
{/i, . . . , fk-i,fk+i,fk, fk+2, ■ • • , f\F\} generates w. First we note that ip{Fk-i) = 
ip{Fk-i) and tjj{Fk+i) — V'(-Ffe+i)- By assumption, wc know Vc ^ Wc, frf ^ Wd and 
vi = wi for all I £ [n]\ {c, d}. Thus 

{v-w), = vc-wc = ^{Fk+i) - ^{Fk) - (^P{Fk) - ij{Fk-i)') 

- ^'(^^fe+i) - ^P{Fk) - i^{Fk) + V(-Ffc-i) 

and 

-^)d^Vd~Wd = i^{Fk) ~ i^{Fk-i) - (ij{Fk+i) ~ ipiPk)) 
= ^{Fk) - ^{Fk-i) - (iAFk+i) - i^iPk)) 
= -4>{Fk+i) + HFk) + ^Fk) - ^(Ffc-i). 

Therefore (v — w)c = — (v — w)^ and w — v is parallel to e^ — ec. Moreover, 
c £ supp(v) since w,v > by assumption that v,vir G 'P(V')- Thus v^r — v, a 
generator of C, is parallel to a vector in i?supp(v)- D 

2.3. The construction of the multivariate rational generating function. 

From the knowledge of triangulations of tangent cones of matroid polytopcs (and 
independence matroid polytopcs) wc will recover the multivariate generating func- 
tions. The following lemma is due to [9] and independently [30]. A proof can also 
be found in [3] and [5]. 

Lemma 20 (Brion-Lawrence's Theorem). Let V be a rational polyhedron and V{V) 
be the set of vertices ofV. Then, 

gviz) = 3Ct,(v)(z), 

vl£V{V) 

where C-p{'v) is the tangent cone of v. 
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Thus, we can write the muhivariate generating function of V by writing aU 
muhivariate generating functions of all the tangent cones of the vertices of V. 
Moreover, the map assigning to a rational polyhedron V its multivariate rational 
generating function gvi^) is a valuation, i.e., a finitely additive measure, so it 
satisfies the equation 

.9-PiuP2(z) = gvi + gv2 - gvinVi^ 

for arbitrary rational polytopes Vi and 7^2, the so-called inclusion-exclusion prin- 
ciple. This allows to break a polyhedron V into pieces Vi and 7^2 and to compute 
the multivariate rational generating functions for the pieces (and their intersection) 
separately in order to get the generating function g-p. More generally, let us denote 
by [V] the indicator function oiV, i.e., the function 



1 if X G 

otherwise. 



Let X^ie/ ^ii'Pi] = be an arbitrary linear identity of indicator functions of rational 
polyhedra (with rational coefhcients e^); the valuation property now implies that it 
carries over to a linear identity X^ie/ 9'Pi (^) = of rational generating functions. 

Now let C be one of the tangent cones of V, and let T be a triangulation of C, 
given by its simplicial cones of maximal dimension. Let T denote the set of all 
(lower-dimensional) intersecting proper faces of the cones Ci G T . Then we can 
assign an integer coefficient to every cone Ci £ T, such that the following identity 
holds: 

This identity immediately carries over to an identity of multivariate rational gen- 
erating functions, 

3c(z) = E + E £»-9c.(z)- (5) 

Remark 21. Notice that formula (5) is of exponential size, even when the tri- 
angulation T only has polynomially many simplicial cones of maximal dimension. 
The reason is that, when the dimension n is allowed to vary, there are exponen- 
tially many intersecting proper faces in the set T . Therefore, we cannot use (5) to 
compute the multivariate rational generating function of C in polynomial time for 
varying dimension. 

To obtain a shorter formula, we use the technique of half-open exact decomposi- 
tions [28], which is a refinement of the method of "irrational" perturbations [6, 27]. 
We use the following result; see also Figure 5 and Figure 6. 
Lemma 22. 



(a) Let 



Ee.[C.] + E^'[C'] (6) 



he a linear identity (with rational coefficients Ei) of indicator functions of 
cones Ci C M", where the cones Ci are full- dimensional for i G Ii and lower- 
dimensional for i £ 12- Let each cone he given as 

= { X e M" : (b*^-, x) < /or i e J, }. (7) 
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(modulo lower-dimensional cones) 



Figure 5. An identity, valid modulo lower-dimensional cones, cor- 
responding to a polyhedral subdivision of a cone 




Figure 6. The technique of half-open exact decomposition. The 
relative location of the vector y (represented by a dot) determines 
which defining inequalities are strict (broken lines) and which are 
weak (solid hnes). 



Let y e R" 6e a vector such that {h*j,y) ^ for all i E IiU h, j G Ji- For 
i (z Ii, we define the "half-open cone" 

= I X e M'^ : (b*^, x) < for j e J, with (b*^-, y) < 0, 

^ ' ' (8) 

(b*^, x) < for j e J, with (b*^-, y) > |. 



Then 



(b) In particular, let 



^£,;[C,]=0. (9) 



ieii 



[c] - E [C'l + E [c^] (10) 

c.er c,er 

be the identity corresponding to a triangulation of the cone C, where T is the 
set of simplicial cones of maximal dimension and T is the set of intersecting 
proper faces. Then there exists a polynomial-time algorithm to construct a 
vector y G Q" such that the above construction yields the identity 

[c] = E (11) 

which describes a partition of C into half-open cones of maximal dimension. 

Proof. Part (a) is a slightly less general form of Theorem 3 in [28] . Part (b) follows 
from the discussion in section 2 of [281. □ 
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Since the cones in a triangulation T of all tangent cones C-p (v) of our polytopes 
are unimodular by Lemma 7, we can efficiently write the multivariate generating 
functions of their half-open counterparts. 

Lemma 23 (Lemma 9 in [28]). Let C C M" be an N -dimensional half-open pointed 
simplicial affine cone with an integral apex w £ Ia^ and the ray description 



C = 



I V + Ajbj : > for j e J< and > for j £ J< | (12) 

where J< U = {1, • • • , N} and hj G Z" \ {0}. We further assume that C is 
unimodular, i.e., the vectors hj form a basis of the lattice (Mbi + ■ • ■ + Mbjv) H Z". 
Then the unique point in the fundamental parallelepiped of the half-open cone C is 

a = v+Vb„ (13) 



and the generating function of C is given by 

gciz)^ ^ f ^ (14) 

Taking all results together, we obtain: 

Corollary 24. Let r be a fixed integer. There exist algorithms that, given 

(a) a matroid M on n elements, presented by an evaluation oracle for its rank 
function tp, which is bounded above by r, or 

(b) an evaluation oracle for an integral polymatroid rank function ip: 2^"! > N, 

which is bounded above by r, 

compute in time polynomial in n vectors a^ G Z", hij G Z" \ {0}, and v.; G Z" for 
i € I (a polynomial- size index set) and j — I, . . . , N , where N < n, such that the 
multivariate generating function of'P{M), 'P^{M) and 'P{4'), respectively, is the 
sum of rational functions 

gp(z) ^ V ^ (15) 

and the k-th dilation of the polytope has the multivariate rational generating func- 
tion 

2;'^,+(fc-i)v, 

Proof. Lemma 20 implies that finding the multivariate generating function of 'P(M) , 
V^{M) or Viyiji) can be reduced to finding the multivariate generating functions of 
their tangent cones. Moreover, V{M), V^{M) and 7'('0) have only polynomially 
in n many vertices as described in subsection 2.1 or Lemma 16. Enumerating their 
vertices can be done in polynomial time by Lemma 5 or 17. 

Given a vertex v of r{M), V'^{M) or Viip), its neighbors can be computed in 
polynomial time by Lemma 5 or 17. The tangent cone C(v) at v is generated by 
elements in Ra, where Ra is an elementary set for some AC [n]. See subsection 
2.1 or Lemma 19. We also proved in Lemma 7 that every triangulation of C(v) 
generated by elements in Ra is unimodal and Lemma 9 states that any triangula- 
tion of C(v) has at most a polynomial in n number of top-dimensional simplices. 
Moreover, a triangulation of the cone C(v) can be computed in polynomial time by 
Lemma 12. Finally, using Lemmas 22 and 23 we can write the polynomial sized 
multivariate generating function of C(v) in polynomial time. Therefore we can 
write the multivariate generating function of V{M), V^{M) or V{ip) in polynomial 
time. □ 
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2.4. Polynomial-time specialization of rational generating functions in 
varying dimension. Wc now compute the Ehrhart polynomial i{V, k) ~ #(fc'P n 
Z") from the multivariate rational generating function guvi'^) of Corollary 24. This 
amounts to the problem of evaluating or specializing a rational generating function 
5fc-p(z), depending on a parameter k, at the point z = 1. This is a pole of each 
of its summands but a regular point (removable singularity) of the function itself. 
From now on we call this the specialization problem. We explain a very general 
procedure to solve it which we hope will allow future applications. 

To this end, let the generating function of a polytope V C R" be given in the 
form 

9v{z) = Y.^. J i (17) 
iei •i^J=i^ ^ 

where £i G {±1}, a^ S Z", and by S Z" \ {0}. Let s = max^g/ Si be the maximum 
number of binomials in the denominators. In general, if s is allowed to grow, more 
poles need to be considered for each summand, so the evaluation will need more 
computational effort. 

In previous literature, the specialization problem has been considered, but not 
in sufficient generality for our purpose. In the original paper by Barvinok [1, 
Lemma 4.3], the dimension n is fixed, and each summand has exactly Si = n 
binomials in the denominator. The same restriction can be found in the survey 
by Barvinok and Pommersheim [3]. In the more general algorithmic theory of 
monomial substitutions developed by Barvinok and Woods [4], Woods [41], there is 
no assumption on the dimension n, but the number s of binomials in the denomina- 
tors is fixed. The same restriction appears in the paper by Verdoolaege and Woods 
[39, Lemma 2.15]. In a recent paper, Barvinok [2, section 5] gives a polynomial-time 
algorithm for the specialization problem for rational functions of the form 

where the dimension n is fixed, the number s of different binomials in each denom- 
inator equals n, but the multiplicity 7^ is varying. 

We will show that the technique from Barvinok [2, section 5] can be implemented 
in a way such that we obtain a polynomial-time algorithm even for the case of a 
general formula (17), when the dimension and the number of binomials are allowed 
to grow. 

Theorem 25 (Polynomial-time specialization), (a) There exists an algorithm for 
computing the specialization of a rational function of the form 

gp(z) = ^s^ (19) 
iei •i-ij=i^ ' 

at its removable singularity 2,-1, which runs in time polynomial in the en- 
coding size of its data Ei G Q, a^ G Z" for i G I and G Z" for i £ I , 
j = 1, . . . , Si, even when the dimension n and the numbers Si of terms in the 
denominators are not fixed. 
(b) In particular, there exists a polynomial-time algorithm that, given data Si G Q, 
aj elT^foriel and G Z" for i (z I , j — 1, . . . , describing a rational 
function in the form (19), computes a vector X G Q" with (A, by ) ^ for all 
i,j and rational weights Wi^i for i G I and I = 0, . . . , s^. Then the number of 
integer points is given by 



#(7'nZ") = ;^£,gu;,.z(A,a,)'. (20) 



iei 1=0 
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(c) Likewise, given a parametric rational function for the dilations of an integral 
polytope V, 



the Ehrhart polynomial i{'P, k) = ^{kV D Z") is given by the explicit formula 
*(^'^)- E fE-<^'-')™E Q)^M(A,a.-v,)'-"A k"\ (22) 

m=0 Vie/ l=m ^ ^ J 

where M ~ min{s, dim'P}. 

Proof of Theorem 1. Corollary 24 and Theorem 25 imply Theorem (Theorem 1) 
directly. □ 

The remainder of this section contains the proof of Theorem 25. We follow [3] 
and recall the definition of Todd polynomials. We will prove that they can be 
efficiently evaluated in rational arithmetic. 

Definition 26. We consider the function 

H{x,iu---,^s) = f[-. 

r=i ^ ~ expj-xCJ 

a function that is analytic in a neighborhood of 0. The m-th (s-variate) Todd 
polynomial is the coefficient of a;™ in the Taylor expansion 



H{x,^i, ...,£,s)= ^ td,„(^i, . . . ,6)a 



771—0 

We remark that, when the numbers s and m are allowed to vary, the Todd 
polynomials have an exponential number of monomials. 

Theorem 27. The Todd polynomial tdm(^i, . . . , ^s) can be evaluated for given 
rational data ^i,...,^^ in time polynomial in s, m, and the encoding length of 

The proof makes use of the following lemma. 

Lemma 28. The function h(x) = x/ {\ — c^^) is a function that is analytic in a 
neighborhood o/O. Its Taylor series about x = is of the form 

oo ^ 

Hx) = bnx" where 6„ = — - — —rrrCn (23) 
^-^ nl (n + 1)! 

n—O ^ ' 

with integer coefficients Cn that have a binary encoding length of 0{n^\ogn). The 
coefficients c„ can be computed from the recursion 

Co = 1 

■A, Ai + lA nl , (24) 

""-EM)^^\,. + J(;7^-TT)i^- Vn = i,2,.... 

Proof. The reciprocal function h^^{x) = (1 — e^^)/x has the Taylor series 

(-1)" 



j=i 

OO 



h ^{x) = y^ Qnx" with a„ = 

^ — ' 71 + 1 ! 
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Using the identity h ^{x)h{x) = (X]^o ^"^") (S^^o ^"^") ~ obtain the 

recursion 



&o = — = 1 

bn = -(aife„_i + a2fo„-2 H h a„6o) for n = 1, 2, 

We prove (23) inductively. Clearly bo ~ cq = 1. For n = 1, 2, . . . , we have 

c„ = n! (n + 1)! 6„ 

= -n! (n + 1)! (ai6„_i + a26„_2 H h a„6o) 

A 1 



(25) 



(.7 + 1)! in~jy.in-j + iy. 

-Vf iv>+i + ^' c . 

^tt (j + l)!(«-j)! (n-j + 1)! 

Thus we obtain the recursion formula (24), which also shows that all c„ are integers. 
A rough estimate shows that 

\cn\ <n(?i + l)!7i! \cn-i\ < ( (n + 1 ) !) ^ | c„„ 1 1 , 
thus \cn\ < {{n + 1)!)^", so c„ has a binary encoding length of 0(n? logn). □ 
Proof of Theorem 27. By definition, we have 

oo s 

77 (x, a , • . • , 6 ) = E (a , • • • , 6 )a;" = n '^(^^^- ) • 

m=0 j=l 

From Lemma 28 we have 

h{xij) = Y,l3,^rix'' + o{x'^) where = ^T^VnT^" (^^^ 

with integers c„ given by the recursion (24). Thus we can evaluate tdm(ai ■ • ■ I'Cs) 
by summing over all the possible compositions rii + • • • + = to of the order to 
from the orders Uj of the factors: 

td™(a,...,a) = Phni ■ ■ ■ (3s,n^ (27) 
(ni,...,ns)eZ^ 
h^s— ^ 

We remark that the length of the above sum is equal to the number of compositions 
of TO into s non-negative parts, 



C'sim) 



m + s — V' 
s-1 

(to + s — 1) (m + s — 2) . . . (to + s — (s — 1)) 
= (s-l)(s-2)...2.1 

which is exponential in s (whenever m > s). Thus we cannot evaluate the for- 
mula (27) efficiently when s is allowed to grow. 

However, we show that we can evaluate tdm(ai ■ • ■ i^s) more efficiently. To 
this end, we multiply up the s truncated Taylor series (26), one factor at a time. 
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truncating after order m. Let us denote 



Hsix) = Hs-i{x) ■ hix^s) = H{x, 6, • • • , 6)- 

Each multiplication can be implemented in O(m^) elementary rational operations. 
We finally show that all numbers appearing in the calculations have polynomial 
encoding size. Let S be the largest binary encoding size of any of the rational 
numbers <^i,...,^s- Then every /?j_„ given by (26) has a binary encoding size 
0(Sn^ log'^ n). Let Hj{x) have the truncated Taylor series X]™=o'^j."^" + o(x'") 
and let Aj denote the largest binary encoding length of any „ for n <m. Then 



Hj+i{x) = ^ aj+i,„.T" + o{x™-) with aj+i,„ = '^aj,iP, 



'j,n-l- 



1=0 



Thus the binary encoding size of aj+i.„ (for n < m) is bounded by ^j+0(Sm'^ log'^ m). 
Thus, after s multiplication steps, the encoding size of the coefficients is bounded 
by 0(sSm^ log"^ m), a polynomial quantity. □ 

Proof of Theorem 25. Parts (a) and (b). We recall the technique of Barvinok [1, 
Lemma 4.3], refined by Barvinok [2, section 5]. 

We first construct a rational vector A S Z" such that {\,hij) ^ for all 
One such construction is to consider the moment curve A(^) ~ (1, ^, . . . , ^"~^) G 
R". For each exponent vector b^j occuring in a denominator of (17), the function 
fij '■ C ^ (-^(Oi by ) is a polynomial function of degree at most n—l. Since ^ 0, 
the function /y is not identically zero. Hence has at most n—l zeros. By 
evaluating all functions fij for i G / and j = l,...,SiatM = (n— l)s|/| + l different 
values for ^, for instance at the integers ^ = 0, . . . , M, we can find one £, — ^ that is 
not a zero of any fij. Clearly this search can be implemented in polynomial time, 
even when the dimension n and the number s of terms in the denominators are not 
fixed. We set A = \{^). 

For T > 0, let us consider the points = e^"^ = (cxpjrAi}, . . . ,exp{TA„}). We 
have 

n 

= ncxp{TA;6.y7} = exp{T (A, by )}; 
1=1 

since (A, by) ^ for all all the denominators 1 — Zr'^ arc nonzero. Hence 
for every t > 0, the point Zj- is a regular point not only of g{z) but also of the 
individual summands of (17). We have 



'^tTi n;Li(i-z^) 

= hm Ve, exp{r(A,a,)} 



o+t^ n;Li(l-cxp{r(A,by)}) 



Hm r-- exp{r (A, a,)} \[ ^—-^-^^-^ 



" '"exp{T(A,a,)}]^— /i(-r(A,by)) 

'^^ is/ j=i ^ ' 
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where H{x, ^i, . . . , S^si) is the function from Definition 26. We will compute the 
limit by computing the constant term of the Laurent expansion of each summand 
about r = 0. Now the function r exp{T(A,ai)} is holomorphic and has the 
Taylor series 

exp{r (A, a,)} ^ ^ a.^r' + o(t"-) where a,,; = "^p", (28) 
/=o 

and H{t, ^i, • ■ • , ^sJ has the Taylor series 

Si 

(r, 6 , . . . , ) = ^ td„, (a , . . . , 6 )t™ + o{t^' ) . 

m=0 

Because of the factor t^"' , which gives rise to a pole of order Si in the summand, 
we can compute the constant term of the Laurent expansion by summing over all 
the possible compositions Si = I + {si — I) of the order s^: 

= E ^' n--^~a b -^ £ ^'^^'-'^^ ^'^^ (29) 

iei •^ii=i ^ ' 1=0 

We use the notation 

, _,,^,tds,_,(-(A,b,.i),...,-(A,b,.s,)) 
li'j,; = (-1 ' „ . , 7T-r ^ for I G / and Z = 0, . . . , s,; 

these rational numbers can be computed in polynomial time using Theorem 27. We 
now obtain the formula of the claim. 

Si 

5(1) = E E '^''^ ('^' • 

iGl 1=0 

Part (c). Applying the same technique to the parametric rational function (21), 
we obtain 

#(fc7'nZ") = .gfep(l) 

Si 

= ^ £i ^ Wi,i (A, ai + {k- l)vi)' , 
iei 1=0 

Si I 



I I ^A, a; - V, 
\ m J 

i£l 1=0 m=0 

Si / } \ ^ 

\l — m 



= E ( E E (^) - 

m=0 Vie/ l=m ^ ' / 

an explicit formula for the Ehrhart polynomial. We remark that, since the Ehrhart 
polynomial is of degree equal to the dimension of V, all coefficients of fc™ for 
TO > dimT' must vanish. Thus we obtain the formula of the claim, where we sum 
only up to min{s,dim7'} instead of s. □ 

3. Algebraic Properties of /i*-vectors and Ehrhart polynomials of 

Matroid Polytopes 

Using the programs cdd+ [20], LattE [11] and LattE macchiato [26] we ex- 
plored patterns for the Ehrhart polynomials of matroid polytopes. Since previous 
authors proposed other invariants of a matroid (e.g., the Tutte polynomials and 
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the invariants of [7, 35]) we wished to know how weU does the Ehrhart polyno- 
mial distinguish non- isomorphic matroids. It is natural to compare it with other 
known invariants. Some straightforward properties are immediately evident. For 
example, the Ehrhart polynomial of a matroid and that of its dual are equal. Also 
the Ehrhart polynomial of a direct sum of matroids is the product of their Ehrhart 
polynomials. 

Table 1 . Coefficients of the Ehrhart polynomials of selected ma- 
troids in [7, 35] 



Ehrhart Polynomial 



Speyerl 


1, 


21 
5 ' 


343 63 91 77 29 
45 ' 8 ' 18 ' 40 ' 90 


Speyer2 


1, 


135 
28 


3691 1511 88 39 529 89 
' 360 ' 120 ' 9 ' 8 ' 360 ' 420 


BJRl 


1, 


109 
30 


23 59 9 9 
' 4 ' 12' 4' 20 


BJR2 


1, 


211 
60 


125 33 43 43 
' 24 ' 8 ' 24 ' 120 


BJR3 


1, 


83 


2783 199 391 247 61 


20 ' 


360 ' 24 ' 72 ' 120' 180 


BJR4 


1, 


25 


353 101 193 23 53 


6 ' 


45 ' 12 ' 36 ' 12 ' 180 



We call the last two matroids in Figure 2 in [35] Speyerl and Speyer2 and the 
matroids of Figure 2 in [7] BJRl, BJR2, BJR3, and BJR4, and hst their Ehrhart 
polynomials in Table 1. We note that BJR3 and BJR4 have the same Tutte poly- 
nomial, yet their Ehrhart polynomials are different. This proves that the Ehrhart 
polynomial cannot be computed using deletion and contractions, like is the case for 
the Tutte polynomial. Examples BJRl and BJR2 show that the Ehrhart polyno- 
mial of a matroid may be helpful on distinguishing non- isomorphic matroids: These 
two matroids arc not isomorphic yet they have the same Tutte polynomials and the 
same quasi-symmetric function studied in [7]. Although they share some properties, 
there does not seem to be an obvious relation to Spcyer's univariate polynomials 
introduced in [35]; examples Speyerl and Speyer2 show they are relatively prime 
with their corresponding Ehrhart polynomials. 

Our experiments included, among others, many examples coming from small 
graphical matroids, random realizable matroids over fields of small positive char- 
acteristic, and the classical examples listed in the Appendix of [31] for which we 
list the results in Tabic 2. For a comprehensive list of all our calculations visit [12]. 
Soon it became evident that all instances verified both parts of our Conjecture 2. 

By far the most comprehensive study we made was for the family of uniform 
matroids. In this case we based our computations on the theory of Veronese algebras 
as developed by M. Katzman in [24]. There, Katzman gives an explicit equation for 
the /i*-vector of uniform matroid polytopcs (again, using the language of Veronese 
algebras). For this family we were able to verify computationally the conjecture is 
true for all uniform matroids up to 75 elements and to prove partial unimodality 
as explained in the introduction. 

Lemma 29. The coefficients of the Ehrhart polynomial of the matroid polytope of 
the uniform matroid U'^'" are positive. 

Proof. We begin with the expression in Corollary 2.2 in [24] which explicitly gives 
the Ehrhart polynomial of 7^ ([/'"'") as 
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Table 2. Coefficients of the Elirhart polynomials and /i*-vectors 
of selected matroids in [31] 





Tl 


T 


/i*-vector 


Ehrhart Polynomial 




6 


3 


1, 10, 20, 10, 1 




107 21 49 7 7 
30 ' 4 ' 12 ' 4 ' 20 




6 


3 


1, 11, 24, 11, 10 


\ 


18 11 992 
5 ' 2 ' 2 ' ' 5 


Qe 


6 


3 


1, 12, 28, 12, 1 


]^ 


109 23 59 9 9 
30 ' 4 ' 12 ' 4 ' 20 




6 


3 


1, 13, 32, 13, 1 


]^ 


11 16 6 1 
3 ' ' 3 ' 2 ' 2 


Rfi 


6 


3 


1, 12, 28, 12, 1 


\ 


109 23 59 9 9 
30 ' 4 ' 12 ' 4 ' 20 




7 


3 


21, 98, 91, 21, 1 


]^ 


21 343 63 91 77 29 
5 ' 45 ' 8 ' IS ' 40 ' 90 


^7 


7 


3 


21, 101, 97, 22, 1 


]^ 


253 2809 33 193 61 121 
60 ' 360 ' i ' 36 ' 30 ' 360 


P7 


7 


3 


21, 104, 103, 23, 1 


\ 


127 479 69 17 257 7 
30 ' 60 ' 8 ' 3 ' 120 ' 20 


AG(3, 2) 


8 


4 


1, 62, 561, 1014, 449, 48, 1 


\ 


209 1981 881 119 95 499 89 
42 ' 180 ' 60 ' 9 ' 12 ' 180 ' 210 


AG' (3, 2) 


8 


4 


1, 62, 562, 1023, 458, 49, 1 


\ 


299 4007 5401 122 2911 1013 77 
60 ' 360 ' 360 ' 9 ' 360 ' 360 ' 180 




8 


4 


1, 62, 563, 1032, 467, 50, 1 




524 1013 1379 125 743 257 136 
105 ' 90 ' 90 ' 9 ' 90 ' 90 ' 315 


^8 


8 


4 


1, 62, 563, 1032, 467, 50, 1 




524 1013 1379 125 743 257 136 
105 ' 90 ' 90 ' 9 ' 90 ' 90 ' 315 


Qs 


8 


4 


1, 62, 564, 1041, 476, 51, 1 




2099 4097 1877 128 337 1043 61 
420 ' 360 ' 120 ' 9 ' 40 ' 360 ' 140 


Ss 


8 


4 


1, 44, 337, 612, 305, 40, 1 




1021 377 475 193 511 65 67 
210 ' 36 ' 36 ' 18 ' 90 ' 36 ' 252 


Vs 


8 


4 


1, 62, 570, 1095, 530, 57, 1 




2117 4367 2107 146 1133 1133 193 
420 ' 360 ' 120 ' 9 ' 120 ' 360 ' 420 


Ts 


8 


4 


1, 62, 564, 1041, 476, 51, 1 


\ 


2099 4097 1877 128 337 1043 61 
420 ' 360 ' 120 ' 9 ' 40 ' 360 ' 140 




8 


4 


1, 62, 569, 1086, 521, 56, 1 




151 2161 3103 143 1669 559 41 
30 ' 180 ' 180 ' 9 ' 180 ' 180 ' 90 


is 


8 


4 


1, 62, 567, 1068, 503, 54, 1 


1, 


527 529 83 137 134 136 47 
105 ' 45 ' 5 ' 9 ' 15 ' 45 ' 105 


J 


8 


4 


1, 44, 339, 630, 323, 42, 1 


1, 


512 193 83 205 361 17 23 
105 ' IS ' 6 ' IS ' 60 ' 9 ' 84 




8 


4 


1, 62, 565, 1050, 485, 52, 1 


1, 


1051 2071 2873 131 1547 529 277 
210 ' 180 ' 180 ' 9 ' ISO ' ISO' 630 




8 


4 


1, 38, 262, 475, 254, 37, 1 


1, 


135 3691 1511 88 39 529 89 
28 ' 360 ' 120 ' 9 ' 8 ' 360 ' 420 




8 


4 


1, 38, 263, 484, 263, 38, 1 


1, 


169 467 581 91 227 68 68 
35 ' 45 ' 45 ' 9 ' 45 ' 45 ' 315 




9 


5 


78, 1116, 3492, 3237, 927, 72, 1 


1, 


307 137141 3223 37807 211 5743 1889 8923 
56 ' 10080 ' 160 ' 1920 ' 16 ' 960 ' 1120' 40320 


AG{2, 3) 


9 


3 


1, 147, 1230, 1885, 714, 63, 1 


1, 


1453 41749 581 34069 927 4541 239 449 
280 ' 3360 ' 32 ' 1920 ' SO ' 960 ' 224 ' 4480 


Pappus 


9 


3 


1, 147, 1230, 1915, 744, 66, 1 


1, 


729 3573 381 1499 243 49 153 57 
140 ' 2S0 ' 20 ' 80 ' 20 ' 10 ' 140 ' 560 


Non-Pappus 


9 


3 


1, 147, 1230, 1925, 754, 67, 1 


1, 


4379 25951 9287 21967 987 2S55 3701 275 
840 ' 2016 ' 480 ' 1152 ' SO ' 576 ' 3360 ' 2688 


Q3(GF(3)*) 


9 


3 


1, 147, 1098, 1638, 632, 59, 1 


1, 


433 3079 4193 5947 167 601 787 149 
84 ' 252 ' 240 ' 360 ' 16 ' 144 ' 840 ' 1680 




9 


3 


1, 147, 1142, 1717, 656, 60, 1 


1, 


723 49 88 24217 1291 625 821 133 
140 ' 4 ' 5 ' 1440 ' 120 ' 144 ' 840 ' 1440 



Letting r = 2, Equation (30) becomes 

(2fc + n - l)(2fc + n - 2) ■ • ■ (2A: + 1) (k + n - 2)(k + n - 3) ■ ■ ■ (k) 

— 77 

{n-iy. ("-!)! 

We next consider the coefficient of /c""^^^ for < p < n — 1, which can be written 
as 

^-^ (n — j-i) ■ ■ ■ (n ~ in) 

l<jl<-<3p<n-l ^ ^ ■'P' 

Y. 



n — 1 (n — ji) ■ ■ ■ (n — jn-i) 

/I n-1 

2n-p-i y / i y _J_ 

\{n — ji) ■ ■ ■ in ~ j-o-i) ^ n — jr. 



n — 1 ^ (n — ■ ■ ■ (n — 7„_i ) 

l<Jl<---<Jp-l<"-2 ^ ^ ■'P ^' 



Hi 
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n-1 



n — j„ n — 1 
ip=l+Jp-i 



= V ( - 2"-p-i V — 

i<ii< - <jp-i<«-2 ^ ^ ^ ' 

(31) 

It is known that the constant in any Ehrhart polynomial is 1 [36], thus we only 
need to show that Equation (31) is positive for < p < n — 2. It is sufficient to 
show that the square-bracketed term of (31), 

1 / 1 \ 
2„_,_i j2 1 + 7 , (32) 

is positive for < p < n — 2 and all 1 < ji < • • • < jp^i < ?i — 2. We can see 
that y;"^\ , , ^— > 1. Moreover, since < p < n - 2 then 2"-P-i > 2 and 
1 + < 2 since ?i > 2, proving the result. □ 

To present our results about the /i*-vector we begin explaining the details with 
the following numbers introduced in [24], which we refer to as the Katzman coeffi- 
cients: 

Definition 30. For any positive integers n and r define the coefficients A^'^ by 

n(r — 1) 

^^'^'■T* = (l + T+--- + T''-i)". 
We also define the vector A"'*^ as ( A'^'' , A'-!'' , . . . ,A 



( 4"''' /I"''' 

\^^0 ' ^1 ' ■ 



i(r-l) 



Looking at the definition of the Katzman coefficient, we see that AJ"" = (p and 

A^'^ — unless j = 0, in which case we have ^q'^ = 1. 

Below we derive some new and useful equalities and prove symmetry and uni- 
modality for the Katzman coefficients. Katzman [24] gave an explicit equations for 
the /i*-vector of uniform matroid polytopes and the coefficients of their Ehrhart 
polynomials, although he did not use the same language. We restate it here for our 
purposes: 

Lemma 31 (See Corollary 2.9 in [24]). Let V{U^'^) he the matroid polytope of the 
uniform matroid of rank r on n elements. Then the h* -polynomial of T'{U^'") is 
given by 



r-l s j 



EEEE(~ir"^"^ 

s=0 j=0 k=0 l>k 

That IS, for ii*(7'(i7"^''))) = {hl,...,h;), 

r-l s j 



T'. (33) 



^r = EEE(-ir"^"'^(:j(y(y^I^:^;;-.)- 



s=0 j=0 fc=0 ^ 

For r = 2 the h* -polynomial of V{U)"''^ is 

i:(;;)r')-„T. (34, 

The following lemma is a direct consequence of Corollary 2.9 in [24]: 

Lemma 32. Let 'P(C/^'") be the matroid polytope of the uniform matroid of rank 
2 on n elements. The h*-vector of V{U'^'"') is unimodal. 
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The rank 2 case is an interesting example already. Although the /i*-vector is 
unimodal, it is not always symmetric. Next we present some useful lemmas, the 
first a combinatorial description of the Katzman coefficients. 

Lemma 33. For i = 0, . . . , n{r — 1) we have 



A^^^= V ( " ) (35) 

Oao + la, + ... + (r-l)a._i=« V«0 , ^1 , " ' ' , 1^ 

ao+aiH ha,,_i=n 

where oq, . . . , a^-i run through non-negative integers. 
Proof. Using the multinomial formula [37] we have 

n(r — 1) 

^ A^'-'-r =(i + T + --- + r-i)" 

i=0 

n 



^ \ ao, ai, • • • , flr-i / 



r— Ija,. 



ao+QlH |-a,.-i=n 

_ / " \2n0ao + lai+-+(r-l)a^_i 

fiQ+aiH |-ar_i=)i \ ' ' ' / 

By grouping the powers of T we get equation (35). □ 

Next we present a generalization of a property of the binomial coefficients. The 
following lemma relates the Katzman coefficients to Katzman coefficients with one 
less element. 

Lemma 34. 

Av= E ^r''^ (36) 

fe=i-r+l 

where we define Ap^^'^ :— when p < or p > (n — l)(r — 1). 
Proof. 

n(r — 1) 
i=0 

= (f + r + --- + T'-i)" = {1 + T + ■ ■ ■ + T'-^y'-\i + T + ■ ■ ■ + r'-^) 

(n-l)(r-l) \ 

^ Ap^'^r (i + T + --- + T'-i) 

(Ti-l)(r-l) (ri-l)(r-l) (n-l)(r-l) 

i=0 4=0 1=0 

(n-l)(j--l) (n-l)(r-l) + l (n-l)(r-l)+r-l 

I— i— 1 i—r—1 

(n-l)(r-l)+r-l 

i=0 

Thus we get equation (36). □ 

The following lemma relates the Katzman coefficients of rank r with those of 
rank r — f . 
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Lemma 35. 



k=0 




or m 



other words 



= E ur'""- 

k+l=i ^ 
0<k<n 



k+l=i 
0<k<n 
0<i<fc(r-2) 



Proof. From Definition 30 

n(r — 1) 

J2 A^^'^T'={l+T 




□ 

Lemma 36. The Katzman coefficients are unimodal and symmetric in the index 
i. That is, the vector (^Aq'^ , A^'^ , . . . ,A^'^^_^^^ is unimodal and symmetric. 

Proof. We first prove symmetry. Considering equation (35) we assume that 

Oflo + lai + • • • + (r — l)ar-i = i and aq + ai + • • • + a^-i = n. 

These two assumptions imply that 

((r - 1) - 0)ao + {{r - 1) - + • • • + ((r - 1) - (r - l))a,_i 
= {r - l)ao + {r - l)ai + ••• + (?'- l)a„ — Ooq - lai -••■-(?'- l)ar-i 
= (r — — z. 

Therefore 



^i'^ = E 



Oao + laiH hC?"— l)^^--!— i 

ao+oiH har-l=n 



n 

Ho, ai, ■ • ■ , flr-l 



(r-l)ao + (r-2)ai+---+0a^_i=i 
ao+aiH |-a,,_i=n 



V ( 1 = ^"'"^ 

^ 0-0, 0.1, - ■ ■ , a,.-!/ ^ ' 



To prove unimodaUty we proceed by induction on n, where r is fixed. First, A. 

' ^(ri-l)(r-l) 



is unimodal. Assume for n — 1 that A" = (^q ^'^ , A^ , ■ ■ ■ ,A" ^'^ 



is unimodal. Using equation (36) and the fact that A"^^'*^ is symmetric we get 
that A"''" is unimodal. To help see this, one can view equation (36) as a sliding 
window over r elements of the vector A"~ , that is, A"'^ is equal to the sum of 
the r elements in a window over the vector A""^''^. As the window slides up the 
vector A""^'*", the sum will increase. When the window is on the center of A" 
symmetry and unimodality of A"~^'' will imply unimodality of A"'' . □ 
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Now we use the explicit equation for the /i*-vector of uniform matroid polytopes 
to prove partial unimodality of rank 3 uniform matroids. First we note that the 
coefficient of T\ hi , in equation (33) is 



s=o j=a k=o 



A'; 



-k}{r-s)- 



Letting the rank r = 3, and using equation (33), we get the /i*-polynomial (which 
is grouped by values of s from (33)), 



E 

i>0 



An.2 I tn~1.2 ^n— 1,2 
— /±2l +^2/ ' "^2(1-1) 



4Tl— 1,1 , — 2,1 



rpl 



Now using that A"'^ ~ (") and A"'^ = SQ{i), where Sj{p) 
get 



1 if P = j 
else 



E 



A'^l^ + n 



n- 1 
21 



n- 1 
21-2 



+ 1 [Soil] -25o{l) + 25i{l) + Soil) -25iil) + S2il)) 



E 



A3, +n 



n — 1 
21 



n- 1 
21-2 



Using properties of the binomial coefhcients, we see that 



n - f 
21 



n — \ 
21-2 



n — 1 

21 
n - 1 
21-1 

n 



n — 1 
21 - I 
n-l 
21-2 



n-l\ /n-l 
21 I ^ \2l-2 



^21 - 1, 

So the /i*-polynomial of rank three uniform matroid polytopes is 



E 

l>0 



3' ^21-1 



62(1) 



(37) 



Using Lemma 36, the coefficient of T'. if 3/ < n, is 



E 



k+p=3l 
0<p<k<n 



n\ / k\ ( n 



'^2(0 



(38) 



' n ' 
.3^ 



3/ - 1 j I 1 



n 

31 - 131 /2\ 



31- [3l/2\\ 
L3//2J ) 

11 



21 - 1 



Next we show that when g is fixed A^' is a polynomial of degree g in the 
indeterminate n, with positive leading coefficient. Assume g < n. Considering 
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Lemma 35 and when g <n, 

k+p=g ^ ^ 
0<p<k<n 

where is a polynomial of degree q with positive leading coefRcient. The highest 
degree polynomial in the sum is (p , a degree g polynomial. Hence A^'^ is a 
polynomial of degree g in the indeterminate n, with positive leading coefficient. If 
g>n, then A^^^ ^ 

since the Katzman coefficients are symmetric by Lemma 

36. 

Proof of Theorem 3 Part (2). Let / be a non-negative integer. From above we see 
that A^'^ is a degree g polynomial in the indeterminate n, with positive leading 
coefRcient. Equation (37) is the ft,*-polynomial of U^-^ , which is a sum of polyno- 
mials in n, the highest degree polynomial being A^l , a polynomial of degree 3L So, 
h* — h*_^ is the difference of a degree 31 and 3(/ — 1) polynomial, hence hi — hi_^ 
is a degree 31 polynomial with positive leading coefficient. For sufficiently large n, 
call it n(/), hi — hl_^ is positive for < Z < I. Hence, the /i*-vector of C/^'" is 
non-decreasing up to the index I for n > n{I). □ 

One might ask if equation (39) has a simpler form. We ran the WZ algorithm 
on our expression, which proved that equation (39) can not be written as a linear 
combination of a fixed number of hypergeometric terms (closed form) [32]. There 
is still the possibility that this expression has a simpler form, though not a closed 
form as described above. 
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